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OPTICAL  FILTERING  TECHNIQUES 


The  purpose  of  this  study  is  to  examine  the  potential  of 
using  optical  information  processing  technology  for  adaptive 
beamforming.  The  work  reported  here  is  a  continuation  of  the  work 
described  in  the  report  R AD C-TR- 8 5  -  8 0.  One  of  the  conclusions 
that  we  reached  in  the  first  part  of  the  program  is  that  adaptive 
broadband  phased  array  processing  is  an  area  where  analog, 
acoustooptic  signal  processing  technology  could  be  fruitfully 
applied.  Our  work  during  the  second  phase  of  this  program  which 
is  reported  here  has  been  concentrated  in  this  area. 

A  schematic  diagram  of  a  broadband  phased  array  is  shown  in 
Fig.  1.  It  consists  of  N  antenna  elements  with  each  element 
having  n  delayed  taps  in  time.  This  allows  both  temporal  and 
spatial  processing  to  be  performed,  a  capabiltity  that  is 
essential  in  the  broadband  case  for  two  reasons.  First,  we 
require  multiple  taps  in  the  broadband  case  simply  in  order  to  do 
beamforming,  i.e.  point  the  antenna  in  the  direction  where  the 
signal  is  coming  from.  In  addition,  the  multiple  taps  allow  us  to 
place  nulls  in  the  antenna  pattern  in  any  direction  and  any 
frequency  at  which  there  is  unwanted  interference.  The  multiple 
taps  provide  additional  degrees  of  freedom  that  increase  the 
total  number  of  jammers  (distributed  in  frequency  and  angle)  that 
can  be  cancelled  to  nN.  Equally  important,  the  multiple  taps 
allow  us  to  optimally  cancel  any  jamming  waveform  within  the 
bandwidth  of  the  signal.  In  an  adaptive  implementation  of  such  an 
array  the  optimum  nN  weights  should  be  continuously  calculated. 


Figure  1.  Broadband,  adaptive  phased  array. 


This  oan  be  accomplished  by  forming  the  nN  X  nN  covariance  matrix 
of  all  the  signals  at  each  tap  and  each  antenna  element, 
and  multiplying  the  Inverse  of  the  covariance  matrix  with  the 
steering  vector  to  obtain  the  optimum  weights.  In  the  previous 
report  we  have  argued  that  at  present  there  is  no  efficient 
method  for  optically  inverting  matrices  with  suitable  accuracy.  A 
digital  implementation  of  such  a  matrix  inversion  would  also  be 
very  difficult  because  of  the  very  large  covariance  matrices 
involved  in  a  broadband  phased  antenna  array.  We  have  thus 
concentrated  on  the  optical  implementation  of  iterative 
algorithms.  In  this  case  the  optimum  weights  are  estimated  by 
multiplying  the  output  of  the  antenna  with  each  delayed  signal  in 
all  the  antenna  elements,  low  pass  filtering  the  result,  and 
subtracting  it  from  the  steering  signal.  Let  s(t)  denote  the 
output  signal  of  the  antenna  and  let  s^t-jA)  be  the  signal  at 
the  ith  antenna  element  and  the  jth  tap.  Then 

s(t)  =  n  w±( j)si(t-j  A) ,  (1) 

1  J 

where  wA(J)  is  the  weight  at  the  jth  tap  of  the  ith  antenna 
element.  The  weights  are  calculated  as  follows: 

t 

wi(J)  =  S  A(J>  -  G  J  s(t')s1(t'-JA)dt',  (2) 

t-T 

where  S A ( J )  is  the  steering  vector  and  G  is  the  feedback  gain. 
When  equations  (1)  and  (2)  are  simultaneously  satisfied,  i.e. 
when  steady  state  has  been  reached,  the  weights  that  are 


calculated  by  Eq.(2)  are  approximately  equal  to  the  optimum 
weights.  Moreover,  if  none  of  the  eigenvalues  of  the  noise 


oovariance  matrix  are  zero  and  the  feedback  gain  is  sufficiently 
low  then  the  Iterative  loop  will  converge  to  this  optimum 
solution . 

The  operation  that  needs  to  be  performed  to  Implement 
equation  (1)  is  a  multichannel  one  dimensional  convolution  (the 
summation  over  the  index  j)  followed  by  a  summation  across  the 
different  channels  to  produce  the  signal  s(t).  Equation  (2)  is 
also  a  multichannel  system  that  correlates  s(t)  with  the  signals 
received  in  each  antenna  element  and  continuosly  produces 
estimates  of  all  the  weights.  The  weights  are  then  subtracted 
from  the  steering  signal  and  fed  back  to  equation  (1)  for 
continuous  iteration.  This  is  shown  in  block  diagram  form  in 
Figure  2,  with  the  symbol  •  denoting  correlation  and  4^  denoting 
convolution.  The  operations  required  for  the  implementation  of 
the  block  diagram  in  figure  2  are  very  well  matched  with  the 
capabilities  of  optical  signal  processing  systems.  One 
dimensional  correlations  and  convolutions  can  be  readily 
Implemented  with  acoustooptic  devices,  the  only  mature  spatial 
light  modulators  that  are  presently  available.  The  multichannel 
capability  can  be  Implemented  by  making  use  of  the  third 
dimension  that  is  available  in  an  optical  system  through  multiple 
transducer  acoustooptic  devices.  Multichannel  acoustooptic 
devices  have  been  developed  relatively  recently  and  they  are  not 
as  well  understood  and  characterized  as  single  channel  devices  . 
However,  devices  with  up  to  32  channels  and  bandwidth  in  excess 
of  100  MHz  are  commercially  available  and  a  larger  number  of 
channels  is  within  the  capabilities  of  current  technology.  Thus, 


we  estimate  that  multichannel  acoustooptic  technology  is  now  or 
will  be  in  the  near  future  sufficiently  well  developed  for  this 
application.  A  careful  examination  of  the  properties  of 
multichannel  AO  devices  ought  to  be  undertaken  however  in  order 
to  estimate  the  expected  performance  of  a  broadband  acoustooptic 
adaptive  array  processor  that  is  Implemented  with  them. 

Having  identified  in  prinoiple  the  general  approach  for 
building  an  adaptive,  broadband  phased  array  processor  using 
multichannel  acoustooptic  technology,  we  now  turn  to  the  choice 
of  specific  architectures.  There  are  two  basic  methods  to:  perform 
co n v o 1 u t io ns / cor r e 1  at ions  with  acoustooptics:  space  or  time 
integration.  There  are  also  several  possible  variations  of  each 
basic  architecture.  In  making  an  appropriate  choice  for  this 
application  we  must  first  take  into  consideration  the  faot  that 
all  the  input  signals  to  the  correlators  and  the  convolvers 
need  to  be  dynamic.  This  fact  narrows  the  choise  considerably. 
Secondly,  the  correlator  and  convolver  blooks  need  to  be 
compatible;  the  output  of  the  convolvers  becomes  the  input  to  the 
correlators  and  vice  versa.  Space  integrating  systems  produce  the 
output  as  a  function  of  time  whereas  time  integrating  systems 
produce  it  as  a  function  of  space.  Thus  space  Integrating  systems 
are  relatively  more  staightforward  to  use  because  we  can  readily 
address  an  acoustooptic  device  with  the  temporal  signal  that  is 
produced  by  them.  The  output  of  a  time  integrating  processor  is 
typically  detected  as  a  function  of  space  on  a  detector  array. 
The  relatively  slow  read-out  of  the  detector  that  converts  the 


detected  signal  to  a  temporal  electronic  signal  makes  the  use  of 


tine  integrating  processors  that  utilize  detector  arrays 
impractical.  There  is  an  alternative  solution  however  which 
involves  the  use  of  optically  addressed  spatial  light  modulators 
which  can  act  as  time  integrating  detectors  and  simultaneously  as 
light  modulators.  This  allows  the  same  device  to  be  used  both  as 
an  output  device  for  one  system  and  an  input  device  for  the 
second.  We  have  investigated  the  possibility  of  using  nonlinear 
optical  crystals  to  accomplish  this  goal.  The  results  of  this 
investigation,  presented  in  the  paper  attached  as  Appendix  A  to 
the  report,  are  very  encouraging  and  we  have  used  this  technique 
in  the  implementation  of  adaptive  filters  as  we  will  discuss 
later. 

Even  with  the  restrictions  mentioned  above  we  are  still  left 
with  considerable  flexibility.  We  can  select  a  purely  space 
integrating  approach  where  both  the  convolution  and  the 
correlation  are  performed  with  space  integrating  processors  or  a 
purely  time  integrating  approach.  Final  1 y  ,  hybrid  approaches  are 
possible  where  time  integration  is  used  for  one  of  the  processors 
and  spatial  integration  for  the  other.  Actually  three  of  these 
four  possible  combinations  have  potential  merit  and  we  have 
explored  them  in  some  detail.  The  only  possibility  that  does  not 
seem  to  be  useful  is  a  fully  time  integrating  approach. 

The  purely  space  integrating  approach  involves  the  use  of 
counter  propagating  multichannel  acoustooptic  devices  to  perform 
the  convolution  whereas  the  correlations  are  performed  with 
copropagating  acoustooptic  devices.  This  system  is  described  in 


detail  in  the  paper  attached  as  Appendix  B  to  this  proposal.  The 


second  approach  we  have  considered  is  the  use  of  time  integrating 
correlations  combined  with  space  integrating  convolvers.  This 
approach  has  been  previously  investigated  by  Rhodes  and  Penn  who 
utilized  optically  addressed  spatial  light  modulators  (e.g.  the 
Hughes  liquid  crystal  light  valve)  as  the  time  integrating 
detector  for  the  correlator  part  of  the  system  and  also  as  a 
spatial  light  modulator  for  the  space  integrating  convolver 
section.  We  have  investigated  the  possibility  of  using 
photoref racti ve  crystals  for  the  same  purpose.  The  results  of 
this  investigation  are  described  in  the  paper  attached  as 
Appendix  C.  The  advantages  of  photoref ract i ve  crystals  are  very 
high  resolution,  excellent  optical  quality  and  contollable 
integration  time.  A  disadvantage  of  the  photoref racti ve  crystals 
that  is  shared  by  the  optical  spatial  light  modulators  is  the 
relatively  long  response  time  (a  hundred  milliseconds  or  more) 
which  translates  to  long  adaptation  time.  The  emergence  of  new 
ohotorefractive  materials,  such  as  GaAs,  with  much  faster 
response  time  than  the  BSO  crystal  used  in  our  experiments 
promises  to  alleviate  this  problem. 

The  third  and  final  possibility  we  explored  is  the  use  of  a 
time  integrating  convolver  and  a  space  integrating  correlator.  We 
will  consider  a  single  channel  version  of  this  system  that 
performs  only  temporal  adaptive  filtering.  This  system  can  be 
extended  to  spatial  and  temporal  filtering  through  the  use  of 
the  third  dimension  and  multichannel  acoustooptic  technology  in  a 


■anner  completely  analogous  to  the  systems  described  earlier.  We 
will  describe  the  operation  of  the  system  for  the  aotive  case 


where  we  are  looking  for  a  known  signal  s(t)  and  jammers  are 
present.  The  optical  system  is  shown  in  Figure  3.  The  input 
signal  x(t)  is  used  to  temporally  modulate  the  laser  diode  in  the 
time  integrating  portion  of  the  system  and  it  is  also  applied  as 
the  input  to  the  acoustooptic  device  in  the  space  integrating 
portion  of  the  system.  The  input  to  ♦’he  acoustooptic  device  in 
the  time  integrating  part  of  the  system  is  the  signal  s(t)  minus 
the  output  signal  from  the  space  integrating  system  amplified  by 
the  gain  G.  The  time  integrating  correlation  forms  as  an  optical 
spatial  modulation  function  on  the  nonlinear  crystal.  This  signal 
is  then  convolved  with  the  input  signal  x(t)  in  the  space 
integrating  part  of  the  system.  The  output  of  the  system  is  the 
signal  e(x/v)  that  forms  on  the  nonlinear  crystal  and  it  can  be 
written  as 


e(x/v)  = 
x»(t)  = 


/ 


J 


x(t)[s(t-x/v)-Gx'(t-x/v)Jdt 
x( t-x/v) e ( x/v) dx. 


(3) 


Taking  the  Fourier  transform  of  the  above  equations  and  solving 
for  the  transform  of  the  of  the  output  signal  E(w),  we  obtain: 

E ( w )  =  X(w)S*(w)/(1+G|N(w) |2)  <U) 

where  |N(w)|2  is  the  power  spectrum  of  the  noise  at  the  input  and 
we  have  neglected  the  effects  of  the  finite  windows  of  the 
acoustooptic  devices  and  the  finite  integration  time  of  the 


nonlinear  crystal.  This  system  is  in  many  ways  similar  to  the 


■V 


system  described  in  appendix  C  with  the  important  distinction 
that  now  the  primary  filter  is  time  integrating  whereas  the  space 
integrating  system  is  in  the  feedback  loop.  This  allows  very  long 
codes  s(t)  to  be  correlated  as  it  would  be  required  in  a  spread 
spectrum  application. 

The  characteristics  of  the  adaptive  filters  we  have  been 
discussing  are  in  general  excellent  in  terms  of  achievable 
bandwidth,  number  of  antenna  elements  and  number  of  taps  per 
element  (i.e.  total  number  of  degrees  of  freedom),  overall  power 
consumption  and  size.  The  adaptation  speed  is  slow  for  the  time 
integrating  approach  with  currently  available  nonlinear  crystals 
or  optically  addressed  spatial  light  modulators,  but  it  is  very 
fast  in  a  space  integrating  implementation.  The  principal 
remaining  question  is  the  maximum  null  depth  that  can  be  achieved 
with  these  analog  implementations.  There  are  two  main  factors 
that  limit  this  performance  criterion:  the  limited  dynamic  range 
of  the  devices  used  and  instabilities  in  the  feedback  loop  caused 
by  phase  distortions  in  the  optical  system.  This  issue  is 
addressed  in  appendix  C,  but  we  will  briefly  summarize  our 
findings  here.  In  the  space  integrating  approach  the  maximum  null 
depth  that  can  be  achieved  in  the  absence  of  phase  distortion  is 
essentially  equal  to  the  dynamic  range  of  the  detectors  used  in 
the  system.  These  detectors  can  have  in  excess  of  80dB  dynamic 
range; hence  If  we  are  willing  to  use  a  sufficiently  powerful 
light  source  to  utilize  fully  the  detector  dynamic  range,  tlen  we 
can  in  principle  achieve  very  deep  nulls.  In  our  experiments 


however  the  phase  errors  in  the  feedback  loop  cause  oscillations 
to  occur  before  the  dynamic  range  limits  are  reached.  These  phase 
errors  occur  for  entirely  different  reasons  in  the  different 
approaches  but  they  have  precisely  the  same  effeot:  If  the 
feedback  gain  is  sufficiently  increased  we  reach  a  condition  of 
positive  feedback  that  can  sustain  oscillation  at  some 
frequencies  within  the  bandwidth.  Thus  the  first  priority  in 
building  a  system  of  this  type  with  very  high  jammer  suppression 
is  to  find  ways  to  have  a  flat  phase  response  over  the  entire 
bandwidth  through  choice  of  the  architecture  and  clever  design 
rather  than  Increasing  the  dynamic  range  of  the  devices  used. 


appendix  a 


A  BIAS  UU  XUU  IBTP0P1TIBQ  0PI1C1L  C01BELAT01 
PSIBG  A  fBPIOPlf BACTITE  C1ISTAL 


Demetri  Psaltls,  Jeffrey  Tu,  end  John  Roeg 
Department  of  Electrical  Engineering 
California  Institute  of  Technology 
Pasadena,  California  9112S 


Abstract 

An  a c ou s t o -o p t i c  time  integrating  correlator  is 
demonstrated  using  a  p  ho  t  o  r  e  f  r  ac  t  i  v  e  crystal  as  the  time 
integrating  detector. 

I.  Introduction 

Tlae  integration1  has  proven  to  be  a  powerful  technique 

in  optical  signal  processing  and  has  been  used  in  a  wide 

variety  of  a  r  c  h  i  t  e  c  t  u  r  e  s  .  A  major  drawback  of  time 

integrating  processors  is  the  buildup  of  bias  in  additon  to 

the  signal.  This  occurs  because  the  pbotogenerated  charge 

that  is  integrated  on  the  detector  is  proportional  to  the 

intensity  of  the  optical  signal  which  makes  it  necessary  to 

represent  bipolar  signals  on  a  bias.  The  effective  system 

dynamic  range  at  the  output  is  given  by  DR'  *  DR 

(SBR / ( 1 +SBR) )  where  DR  is  the  dynamic  range  of  the  output 

detector  and  SBR  is  the  signal  to  bias  ratio  on  the 
2 

detector.  In  most  cases  of  Interest,  the  SBR  is  much 
smaller  than  unity  and  thus  the  added  bias  significantly 


educes  the  usable 


The  Boat  frequently  need  aethod  for  separating  the 
elgnel  from  tbe  blea  loTolite  placing  the  algnal  on  a  apatlal 
oarrler  and  than  electronical ly  filtering  tbe  output  of  tbe 
Integrator.  Tbla  aethod  of  blaa  reaoTil,  however,  does  not 
solve  tbe  dynaalc  range  problea  aince  tbe  proceaalng  la  done 
after  tbe  detection  of  tbe  algnal.  Also,  an  additional 
constraint  la  placed  on  tbe  revolution  of  tbe  detector,  aince 
tbe  pixel  aeparatlon  auat  be  leas  than  one  half  of  tbe  period 
of  the  carrier  being  recorded,  which  will  reault  In  a 
significant  reduction  In  the  available  apace-bandwldtb  product 
at  tbe  output. 

In  this  paper,  a  new  aetbod  for  performing  time 
Integrating  correlation  la  deacribed  using  a  pbotoref rac ti v e 
Bismuth  Silicon  Oxide  (BSO)  crystal  as  tbe  time  Integrating 
element.  Tbe  correlation  Is  formed  on  a  spatial  carrier  In 
tbe  crystal  and  read  out  with  an  auxiliary  beam.  Since  only 
tbe  signal  recorded  on  a  apatlal  carrier  is  stored  In  tbe 
p bo t ore f r ac t 1 v e  crystal,  tbe  diffracted  light  that  Is 
detected  contains  tbe  correlation  Information  without  tbe 
bias.  Tbe  bias  does  not  reduce  tbe  dynamic  range  of  tbe 
output  detector  used  for  final  readout,  but  rather  tbe 
diffraction  efficiency  of  tbe  BSO  crystal.  In  addition,  tbe 
resolution  of  tbe  BSO  crystal  is  very  much  higher  than  that 
of  a  CCD  detector,  allowing  tbe  correlation  of  very  high 
space-bandwidth  signals  to  be  formed  on  a  carrier.  Finally, 
since  tbe  reault  of  the  tlae  Integrating  correlator  Is  read- 


out  optically,  the  output  can  be  easily  Interfaced  with 
other  optical  systems,  thus  asking  new  architectural  designs 
possible. 

In  section  II.  the  theory  of  optical  recording  In 
photorefractlve  crystal  is  reviewed  and  extended  to  the  use 
of  pbotorefractlve  crystals  as  tlae  Integrating  elements. 
The  architecture  and  experlaental  results  are  described  In 
section  III.  Dynaalc  range,  linearity,  aysten  limitations, 
and  other  performance  aspects  are  discussed  In  section  IV. 

II.  Photorefractlve  Crystals  as  Tlae  Integrating 
Optical  Detectors 

When  a  photorefractlve  BSO  crystal  is  Illuminated  by  an 
intensity  grating,  electrons  are  excited  from  traps  Into  the 
conduction  band.  These  charges  migrate  due  to  diffusion  and 
drift  from  an  externally  applied  electric  field  and  then 
recombine  in  dark  regions,  creating  a  spatially  varying 
Internal  space~charge  field.  This  field  modifies  the  index 
of  refraction  in  the  crystal  through  the  linear  electrooptic 
effect  and  as  aresult,  a  holographic  phase  grating  is 
recorded  In  the  crystal.  Grating  formation  in 
photorefractlve  media  has  been  extensively  studied  and 
modeled**4.  We  will  show  here  that  the  pbotorefractlve 
orystal  acts  as  a  time  Integrating  element. 


Let  the  Intensity  incident  on  the  crystal  be  as  follows 


I(x.t) 


I0  ♦  lidjU.t)*111)  for  t>0 
0  otherwise 


(1) 


iBiualDg  that  aalf  diffraction  effects  are  negligible  and 
that  tbe  spatial  variations  of  I^fx.t)  are  small  compared  to 
tbe  grating  frequency  k.  then  tbe  Intensity  of  tbe  llgnt  that 
la  diffracted  when  tbe  crystal  la  illualnated  by  a  read-out 
bean  can  be  abovn  to  be:3 


I0ut(x't) 


i--i-  .-»/«  fSiiizi 

lx  J0  Jn 


dt 


I 

I  I 
I 


Ip  la  tbe  intensity  of  tbe  readout  beta  and  la  a  coaplex 
conatent  Involving  tbe  aaterial  paraaeters  of  tbe  crystal, 
tbe  grating  frequency,  and  tbe  applied  electric  field.  x  la 
tbe  complex  time  constant  of  tbe  apace  charge  field  and  la 
given  by3  t  *  I2/I0-  I2  18  *l0°  *  coaplex  constant  that 
depends  on  tbe  p b o t o r e f r a c t 1 v e  aaterial  used  and  tbe 
experimental  conditions.  Iq  la  tbe  average  light  Intensity 
Incident  on  the  crystal  during  exposure. 

If  Ij (x, t)  la  expanded  Into  Its  temporal  Fourier 
components . 


r- 

Ij(x.t')  -  jljCx.w) 


eiwt#  dw. 


then  tbe  output  Intensity  for  t>>x  can  be  written  as  follows: 


ou  t 


(x.t)  =  I 


I  K 


Ij (X . w) 


I  xlfi  J- 


,  1  w  t 


dw  I  I, 


(1/x  +iw) 


frequency  1/|t|,  which  la  appro* lma tel y  equivalent  to  the 
output  of  a  eliding  window  Integrator,  with  Integration  tlae 
x.  Tbua, 


i  *1  f 

Iout<x*t>  S  • - J 

I  X  J 


I  K.  ft+x  I,(x,t’>  I 


dt'  I  IR 

i  * 


Hence,  the  output  Intenalty  can  be  treated  as  the  square  of 


the  normalized  integration  of  tbe  algnal  Ij.  An  interesting 
observation  la  that  if  Ij(x,t)  were  Independent  of  tine,  then 
tbe  output  intensity  depends  only  on  tbe  ratio  of  tbe 
aodulated  Intenalty  Ij  to  tbe  DC  intensity  Ip. 

III.  Experimental  Demonstration 

|  A  schematic  diagram  of  tbe  experimental  system  Is  shown 

in  Pig  1.  Tbe  input  electrical  signals  are  mixed  with  the 
center  frequency  of  tbe  acousto-optic  devices  (AODs)  and  fed 
I  into  the  AODs.  Tbe  first  AOD  is  illuminated  by  a  collimated 

wave  and  the  upahlfted  diffracted  order  is  Imaged  onto  tbe 
aecond  AOD,  then  relmaged  onto  the  p h o t or ef rac t 1 v e  BSO 
I  crystal.  The  aecond  AOD  is  oriented  such  that  tbe  acoustic 

signal  is  counterpropagating  with  respect  to  tbe  image  of  the 
acoustic  signal  from  the  first  AOD.  The  undiffracted  light 
transmitted  through  the  first  AOD  Is  incident  at  tbe  Bragg 
angle  of  the  aecond  AOD.  The  upsblfted  diffracted  order  of 
»  the  second  AOD  is  also  imaged  on  to  the  BSO  crystal.  The 

a 

I  undiffracted  light  la  spatially  filtered  before  reaching  tbe 


BSO  crystal.  In  this  arrangement,  the  AODs  are  parallel  to 
each  other  .  but  the  diffracted  orders  propagate  at  an  angle 
with  respect  to  each  other  even  though  both  diffracted  beans 
are  temporally  upahlfted.  This  causes  the  signals  from  the 
two  AODs  to  Interf erometrlca 1 1 y  record  the  correlation 
signal  on  the  BSO  crystal  at  a  high  spatial  frequency.  Let 
the  Inputs  to  the  AODs  be 

TjCt)  -  a(t)elwOt 
▼2(t)  -  b(t)elwOt. 

where  wq/2r  Is  the  center  frequency  of  the  AOD.  The 
Intensity  incident  on  the  photorefractlve  crystal  Is 

I(x.t)  ■  | a ( t-x/v) e*TX  ♦  b( t+x/v ) e-lT* | 2 
-  I  a ( t-x / v ) | 2  ♦ | b(t+x/v) | 2  + 

2  fie l a ( t-x / v) b* ( t+x / ▼ ) ei2^x ) ,  (5) 

where  v  Is  the  acoustic  velocity  of  the  AOD  and  y  ■  w 0/v.  We 
will  treat  the  case  where  |a(t)|2  and  |b(t)|2  can  both  be 
approximated  as  constants,  as  la  the  case  for  FM  signals. 
Then,  the  Intensity  pattern  can  be  separated  into  a  DC  term 
Iq  and  a  signal  term  Ij(x,t)  modulating  a  spatial  carrier 
oos(2yx)  in  the  form  of  Eq.  1.  This  intensity  pattern 
results  In  the  formation  of  a  hologram  on  the  photorefractlve 
crystal  as  described  in  the  previous  section.  The  hologram 
la  readout  with  an  auxilary  beam  and  Is  imaged  onto  a  charge 
coupled  device  (CCD)  detector  for  read-out. 

If  the  assumption  that  Ij(x,t)  has  spatial  frequencies 


which  arc  snail  ooapared  to  tha  oarriar  fra quency  2y  la  valid 
than  va  can  usa  Bq.  4  to  obtain  an  azpraaaion  for  tha  output 
intensity  datactad  by  tha  CCD 


i  i,  i2  ift*i 

2o»t<*>  S  2»  I—'!  Ij 

I  v  I  r  t 

and  by  defining  variable  t j 


a(t-x/v)  b*(t+x/v) 

- dt 

I  a  I  2  ♦  Ibl2 

•  t-x/v 


I2 

I. 
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Iout<*> 


I ft+T-X/V 

Ij  aUjlb  (t^x/v)  dt, 

I ' t+X / V 


I2 

I 
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(6) 


Hence  the  ayatei  produces  the  aagnitude  square  of  the 
correlation  between  the  signals  a ( t )  and  b(t)  Integrated  over 
a  finite  interval  x. 

Flint  glass  acousto-optic  cells  driven  at  a  center 
frequency  of  70MHz  were  used  in  the  experiment.  A  symmetric 
linear  chirp  signal  with  bandwidth  Af  •  5Mhz  was  fed  into 
each  cell  to  produce  the  autocorrelation  peak.  The  Bragg 
angle  of  the  AODs  was  0.2°.  which  corresponded  to  a  grating 
frequency  equal  to  35  lines  per  all  llaeter  in  the  BSO 
crystal . 

The  BSO  crystal  used  in  the  experiment  was  cut  in  the 
< 1 1 0  >  direction  and  measured  15  x  15  x  2  am.  An  external 
electric  field  of  7  kV/ca  was  applied  in  the  <001>  direction 
of  the  crystal  which  was  also  the  direction  of  the  grating 
vector . 


The  correlation  was  recorded  on  the  crystal  with  an 
Argon  laser  at  a  wavelength  of  514  na  with  average  intensity 


•qual  to  lpW/ca  .  The  correlation  was  read-out  with  a  Be-Ne 
laser  (X  -  633  nm)  with  lntenalty  ISO  pW/em2.  Cylindrical 
lenaea  (not  shown  In  Flg.l)  were  uaed  to  expand  the  output  of 
the  AODa  thereby  Illuminating  the  full  aperture  of  the  BSO 
cryatal  and  also  to  focua  the  diffracted  light  onto  a  one 
dlnenalonal  CCD. 

The  output  signal  to  bias  ratio  of  a  conventional  tine 
Integrating  correlator  la  reduced  when  the  levels  of  the  two 
signals  are  unequal  and/or  If  there  la  additive  noise 
present  In  the  system.  Both  conditions  were  simulated 
experimentally.  Noise  was  simulated  by  adding  a  70MHz  signal 
to  the  Input  of  one  of  the  AODs.  The  output  of  a  standard 
time  Integrating  correlator  (l.e.  the  correlation  formed 
directly  on  the  CCD),  for  the  noise  free  case  and  equal 
amplitude  signals  Is  shown  In  Fig  2.  This  condition  provides 
the  maximum  signal  to  bias  ratio  for  the  system.  Ve  can  see 
In  Fig  2  that,  there  Is  still  a  strong  bias  term  added  to  the 
correlation  peak.  The  correlation  produced  by  temporally 
Integrating  on  the  photorefractl ve  crystal  Is  shown  In  Fig. 3. 
In  this  case,  all  the  bias  due  to  temporal  Integration  Is 
removed,  and  any  residual  bias  Is  due  entirely  to  dark 
current  from  the  CCD.  The  outputs  of  the  bias  removal 
correlator  with  Input  signal  to  noise  ratios  of  OdB  and  -lOdB 
are  abown  In  Figs  4  and  3,  respectively.  Again,  bias  levels 
which  appear  in  the  figures  were  entirely  due  to  the 


integration  of  dark  current  In  the  output  detector. 


In 


practice,  tba  detector  dark  current  can  be  minimized  by 
increasing  the  intensity  of  the  readout  baas,  thereby 
decreasing  the  required  integration  time  of  the  output  CCD 
detector  and/or  cooling  the  detector. 

III.  Performance 

The  experimental  results  described  In  the  previous 
section  show  a  dramatic  qualitative  improvement  in  the 
correlation  that  la  obtained  when  the  photoref rac ti ve  crystal 
Is  used  Instead  of  the  CCD.  In  this  section,  we  examine 
certain  characteristics  of  this  method  which  are  useful  for 
quantitatively  evaluating  its  performance.  Specifically,  we 
examine  the  linearity.  Integration  tlae,  dynamic  range,  and 
sensitivity  of  the  correlator. 

1)  Linearity 

In  a  conventional  time  Integrating  correlator  (coherent 
or  Incoherent),  the  output  correlation  Is  basically 
proportional  to  the  signals  applied  to  the  AODs. 
Non  1 lnear 1  ties  occur  only  when  we  exceed  the  linear  dynamic 
range  of  the  devices  used,  i.e.  if  the  diffraction  efficiency 
of  the  AOD  exceeds  several  percent  or  the  Integrating 
detector  Is  driven  to  saturation.  In  the  pho toref rac ti v e  time 
Integrating  processor,  the  output  Intensity  is  a  nonlinear, 
monotonies  1 1 y  Increasing  function  of  the  Input  voltage.  The 
nonlinearity  arises  because  of  the  square  law  detection  at 


the  final  rud-out  stage  and  the  recording  mechanism  In  the 
photorafraoti ve  crystal.  The  nonlinear  relationship  la  now 
studied  analytically  and  experlaental  verification  of  the 
theoretical  results  is  presented. 

Let  v^(t)  *  s(t)  be  a  fixed  reference  signal  and  Vj(t)  * 
as(t)  be  an  input  signal  of  varying  aaplltude  (0<a<l).  Since 
the  correlation  term  contains  spatial  frequencies  which  are 
auch  lower  than  the  grating  frequency,  near  the  correlation 
peak  (x«0)  the  intensity  incident  on  the  ph o tor e free t 1 v e 
crystal  la 


l(x.t)  -  (1  ♦  a2  ♦  2a  coa  kx )  |e(t)|2. 


Dalng  Eq.  1.  the  output  intensity  at  the  CCD  la  proportional  to 


Lout 


I  2s  I2 

I - r-  *  I 

I  1  ♦  a2  I 


The  modulation  depth  of  the  intensity  incident  on  the  BSO 
crystal  la 


2a 


a  • 


1  ♦  a' 


and  hence 


‘out 


i2.  -  4a2/  ( 1+a2)  . 


(7) 


fig  6  la  a  graph  of  the  output  Intensity  at  the  correlation 


Bilk  w  a  .  tk*  mdn  1  i t <  an  riant h  tnntriant  nn  Innlriant  nn  tha 


crystal.  The  experlaental  result  la  in  excellent  agreement 
with  the  square  law  relationship  predicted  by  Eq.  7. 

A  plot  of  the  output  intensity  as  a  function  of  the 
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amplitude  of  the  input  signal  is  shown  in  Pig  7.  The 
nonlinear  relationship  between  the  input  and  output  signals 
is  generally  s  disadvantage  since  the  scaling  of  signals  of 
varying  amplitudes  will  be  nonlinear.  This,  however,  will 
not  cause  a  problem  if  the  correlator  is  used  only  as  a 
signal  detection  device,  since  correlation  peaks  will  still 
be  discernible  and  only  the  threshold  level  need  be  adjusted 
accordingly  to  maximize  the  probability  of  detection. 


11)  Integration  Tine 

In  a  conventional  time  integrating  correlator,  the 
integration  tlue  is  Halted  by  the  dark  current  build-up  on 
the  output  detector,  typically  up  to  several  hundred 
ni 1 1 iseconds.  When  the  pbotoref rac t 1 v e  crystal  is  used,  the 
Integration  tine  is  determined  by  the  rise  time  of  the 
Internal  space  charge  field  which  can  easily  be  made  much 


IM 


t  m 


longer.  The  correlation  can  be  read-out  at  any  rate  that  is 


convenient  by  an  auxiliary  detector  array. 

The  Integration  time  Is  approximately  equal  to  |x|. 


where 

I  K,  I 

M  -I— --I  (8) 

I  I0  I  • 

Hence,  the  integration  time  of  the  bias  removal 
correlator  can  be  controlled  by  varying  the  writing 
intensity.  This  control  is  Important  since  the  Integration 


time  can  be  matched  to  the  length  of  the  reference  signal 


V. 


thereby  increasing  the  probability  of  detection  of  a  weak 
signal . 

The  tiae  reponae  of  the  correlation  peak  for  different 
values  of  average  incident  intensity  is  shown  in  Fig.  8. 
Fig.  9  is  a  plot  of  intensity  vs  the  inverse  of  the 
experliental ly  observed  rise  tiae.  There  is  excellent 
agreement  between  the  experlaent  and  equation  8. 

The  integration  tiae,  however,  has  a  finite  range  over 
which  it  oan  be  adjusted.  The  aaxlauB  integration  tiae  is 
Halted  by  the  theraal  effects  in  the  crystal.  If  the  rate 
at  which  carriers  are  generated  theraally  becoae  coaprable  to 
the  rate  at  which  they  are  p h o t o -g e n e r a t e d ,  then  the 
aodulatlon  depth  of  trap  density  will  be  reduced.  As  a 
result,  the  diffraction  efficiency  of  the  grating  will 
decrease.  In  practice,  the  alniaua  integration  tiae  is 
Halted  by  the  aaxlaua  light  intensity  that  is  available  for 
recording.  The  integration  tiae  can  be  reduced  to  30ns,  if 
the  incident  intensity  is  aade  equal  to  18mW/ca  .  This  power 
level.  however.  is  slaply  not  practical  for  aost 
applications . 

ill)  Dynaalc  Range  and  Sensitivity 

Since  the  output  of  the  bias  removal  correlator  is 
presented  without  bias,  the  output  dynaalc  range  of  the 
systea  is  essentially  equal  to  the  dynaalc  range  of  the 
readout  detector  array.  In  order  to  characterise  the 
performance  of  the  systea  we  need  to  deteralne  bow  the  input 


signal  levels  are  sapped  to  this  output  dynamic  range.  Let 
the  dynaalc  range  of  the  photorefractlve  crystal  be  defined 
as  ORgsO  “  ■Bax/»B|n,  where  >Baz  Is  the  aazlaui  modulation 
depth  (m^^-1),  8nd  bb1q  18  the  minimum  Modulation  depth  for 
which  a  diffracted  signal  Is  detectable  above  the  output 
scatter  and  noise  level  of  the  aystea. 

Given  two  input  signals  Vj(t)  ■  a  s(t)  and  v2(t)  *  s(t), 
the  aodulatlon  depth  of  the  light  Incident  on  the  crystal  Is 
a  -  2a/(a2+l).  Thus,  the  alnlaua  detectable  Input  signal  is 
given  by 

amin  "  “min/2  “  1/DRBS0 

The  useful  range  over  which  'a'  can  vary  Is  Halted  by  DRbsq. 
From  aBin,  one  can  define  an  input  dynamic  range  given  by 

DRinput  *  2^aain  *  *^®mln 

The  aost  important  parameter  in  determining  the  system  dynamic 
range  is  the  alnlmum  detectable  aodulatlon  depth  Bm^0> 
Experimentally,  we  measured  the  dynamic  range  to  be  equal  to 
23dB.  This  corresponds  to  a  minimum  aodulatlon  depth  of 
0.142.  We  expect  that  through  careful  design  this  can  be 
substantially  laproved.  However,  all  the  mechanisms  that 
determine  mBln  are  not  fully  understood.  It  Is  believed  that 
besides  detector  noise  and  scattering  froa  the  crystal,  the 
aodulatlon  depth  Is  limited  by  thermal  effects  In  the 
material  and  shot  noise  arising  froa  the  Internal  currents. 

Another  laportant  aspect  of  the  correlator  system  Is  its 
sens  it iv i t v  or  the  alnlaua  signal  to  noise  ratio  that  Is 
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detectable.  This  parameter  la  also  determined  by  tbe  minimum 
detectable  modulation  depth,  a>in  .  Given  a  reference  algnal 
Vj(t)  -  a  a(t)  and  an  Input  algnal  contaalnated  by  additive 
noise,  v2(t)-  b  a(t)  ♦  n(t),  the  modulation  depth  of  the 
intensity  Incident  on  the  crystal  Is 
2ab  |s(t)|2 
<a2+b2) |s(t) |2+cn2  . 

Tbe  reference  level  which  maximizes  a  Is  given  by 
a  -  (b2+en2/|s(t) l2)1/2, 
corresponding  to  a  modulation  depth  of 

b 

<b2+on2/|a(t) I2)2'2. 

In  practice,  optlalzlng  the  reference  level  can  easily  be 
achieved  by  setting  tbe  power  of  tbe  reference  equal  to  the 
total  average  power  of  the  Input  signal. 

Normalizing  the  algnal  and  noise  teras  such  that  I s ( t )  1 2 
■  on2  -1  we  obtain: 

b 

(b2+l)2'2  . 

Thus,  the  alnlaua  Input  SNR  that  produces  a  detectable 
correlation  peak  at  the  output  Is 

r  *.m2- 

Proa  tbe  experiments  1 1 y  measured  value  of  i|io,  tbe 
correlator  should  have  had  a  sensitivity  of  -17dB.  Rowever, 
experlaental  results  showed  a  sensitivity  of  -14dB. 


IT.  Comclaaiom 


Tbe  pho toref rac t 1 v e  time  Integrating  prooessor  that  bas 
been  described  bas  several  advantageous  features:  bias 
removal.  increase  in  tbe  output  space  bandwidth  product,  and 
tbe  ability  to  directly  interface  tbe  result  of  tbe  tine 
integrating  processor  with  other  optical  systems.  Bias  free 
correlation  is  desirable  because  it  allows  us  to  Increase  tbe 
dynamic  range  and  hence  the  sensitivity  of  time  integrating 
processors.  In  tbe  implementation  described  in  this  paper, 
however,  tbe  square  law  detection  at  tbe  output  reduces  tbe 
available  overall  dynamic  range.  A  definite  Improvement  in 
dynamic  range  can  be  obtained  if  tbe  correlation  that  Is 
formed  in  the  pbotorefracti ve  crystal  is  interf erometrica 1 ly 
detected  on  tbe  output  detector.  Another  limitation  of  tbe 
system  described  here  is  tbe  long  integration  time  (several 
seconds).  In  some  applications  this  long  integration  time  is 
desirable  and  could  result  in  extremely  good  sensitivity 
(detection  of  signals  with  very  low  SNR).  However  it  is 
certainly  desirable  to  be  able  to  decrease  tbe  integration 
time  to  several  milliseconds.  This  could  be  accomplished  by 
increasing  tbe  optical  power  of  tbe  writing  beams,  but  this 
is  in  general  an  impractical  solution.  Another  limitation  of 
this  technique  is  tbe  relatively  low  diffraction  efficiency 
that  can  be  obtained  with  BSO  crystals  (2-3%),  wblcb 
reduces  tbe  overall  light  efficiency.  Materials  with  higher 
electrooptic  coefficients,  such  as  Barium  Tltanate.  can 


provide  better  efficiency^  however,  the  t lee  oonstent  that  is 
obtained  with  this  particular  aaterlal  Is  auch  longer  than 
BSO.  New  photorefractl ve  aaterlals  that  are  ourrently  being 
developed  and  show  proalse  for  a  large  laproveaent  in  optical 
sensitivity  as  well  as  higher  electrooptic  coefficients  aay 
provide  a  very  substantial  laproveaent  In  perforaance,  and 
specifically,  reduce  the  total  optical  power  that  is 
required . 
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Figure  2.  Output  of  a  Standard  Time  Integrating 
Correlator  without  Noise 


Figure  3.  Output  of  the  Bias  Removal  Correlator 
without  Noise 
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Ab&lL&Qt 

Space-integrating,  acoustooptic  processors  for  adaptive,  temporal  filtering  are 
examined.  The  basic  architecture  is  then  extended  to  the  space-time  domain  for  application 
in  broadband  phased  array  processing.  An  acou6tooptic  processor  capable  of  such  2- 
dimensional ,  adaptive  processing  is  described. 

A  major  portion  of  optimum  filtering  theory  concerns  itself  with  the  efficient 
separation  of  useful  signals  from  additive  noise.  Fixed  optimum  filters,  such  as  the 
Wiener  filter,  are  applicable  when  the  signal  and  noise  statistics  are  stationary  and  known 
apriori;  the  lack  of  such  apriori  knowledge,  however,  motivates  the  implementation  of 
adaptive  filters  which  estimate  the  required  signal  and  noise  characteristics.  The 
implementation  of  such  filtering  techniques  requires  a  processor  which  must  be  able  to 
compute  various  correlation  functions  of  signal  and  noise  and  change  its  transfer  function 
accordingly . 

The  situation  becomes  considerably  more  complex  when  one  is  required  to  adaptively 
filter  signals  in  the  space-time  domain,  as  in  the  case  of  broadband  array  processing.  The 
number  of  broadband  jammers  that  an  adaptive  processor  can  cancel  without  compromising  its 
directional  discrimination  can  be  used  as  its  perform. ance  m.easure.  Electronic 
implementations  of  such  array  processors1  have  exhibited  limited  performance  with  respect 
to  thi6  measure.  The  transversal  filters  responsible  for  temporally  filtering  the  outputs 
from  the  array  elements  are  of  low  order  and  hence,  operation  is  usually  limited  to 
narrowband  signals,  due  to  hardware  limitations. 

The  bandwidth  requirements  and  the  parallel  nature  of  array  processing  make  optical 
implementations  attractive.  Various  optical  im.pl  ementat  ions  have  been  explored  by 
researchers  ir.  the  area3,3.  In  a  recent  paper4,  we  described  two  optical  adaptive  filters 
for  use  ir,  the  time  domain.  These  implementations  are  strictly  one-dimensional,  requiring 
only  one  dimensional  devices,  permitting  an  extension  to  the  space-time  domain  through  the 
use  of  mu  1 1 1  -  char,  ne  1  AOD's.  In  this  paper,  we  briefly  review  the  operation  of  the 
temporally  adaptive  acoustooptic  filter.  After  this  foundation  has  beer,  established,  an 
adaptive  array  processor  utilizing  multi-channel  AOD's  in  a  space  integrating  architecture 
will  be  described. 

U-i  Adapts  IfiiEpcxAl  Eilttis 

To  be  adaptive,  a  processor  that  optimizes  either  the  mean-square  error  or  SNF  criterion 
must  compute  various  correlation  functions  of  the  input  signal  and  noise  and  vary  its 
filter  response  characteristics  accordingly.  Specifically,  we  will  consider  situations 
where  the  6ignal  is  broadband  and  the  additive  noise  consists  of  strong,  narrowband  jammers 
whose  frequencies  are  unknown;  the  jammer  frequencies  must  be  estimated  for  effective  noise 
rejection . 

Shown  in  Fig. 3  is  a  system  diagram  of  the  Passive  Processor  which  was  shown  to 
adaptively  perform  an  approximate  Wiener  filtering  operation4.  The  operation  of  the  syster 
car.  be  explained  heu r  1  st ica 1 1 y  in  the  following  manner.  Suppose  that  the  input  consists 
aoley  of  a  single  6inusoid,  and  consider  the  feedback  signal,  x(t).  with  the  loop  opened  at 
the  summing  junction.  The  input  sinusoid  is  first  convolved  with  a  replica  of  itself, 
delayed  by  e  to  produce  another  sinusoid  of  the  same  frequency,  also  delayed  by  o,  at  the 
input  of  the  correlator.  This  is  then  correlated  with  the  input  delayed  by  c  to  produce  a 
ainuacid  with  the  same  phase  as  the  input  sinusoid  appearing  at  the  summing  junction,  the 
delays.  o.  cancel  out  due  to  the  cascaded  cor r e 1  at l on  -  con v o 1 u t i on  operations.  The 
possibility  of  statle  cancellation  is  now  apparent  since  the  two  signals,  the  input  and 
feedback  signals,  which  are  subtracted  at  the  junction,  are  identical  in  frequency  and 
phase.  Cancellation  does  not  occur  for  a  broadband  input,  because  it  correlates  poorly  wit: 
delayed  versions  of  itself,  resulting  in  a  negligible  feedback  signal.  £(t>.  The  passive 
processor  thus  discriminates  against  narrowband  components  of  the  input  signal,  while 
preserving  the  desired,  broadband  components  of  the  input,  making  it  suitatle  fci  signal 
estimat ion . 
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Pig. 2  Active  Processor 


For  signal  detection,  where  a  known  broadband  signal  is  to  be  detected  airidst  an 
additive  collection  of  narrowband  jammers  of  unknown  frequencies,  a  different  system  is 
required.  Shown  in  Fig.  2  is  a  system  diagram  of  such  a  system,  the  Active  Processor,  along 
with  a  brief,  mathematical  description  of  its  operation;  a  complete  description  can  be 
found  in  reference  a.  As  with  the  Passive  Processor,  the  basic  operation  of  this  system  can 
be  understood  using  heuristic  arguments.  Suppose  that  the  input  consists  soley  of  a  large 
amplitude  sinusoid  of  frequency  fp.  and  the  desired  waveform.  s(t).  is  a  broadband  signal. 
Kith  the  loop  opened  at  the  summing  junction,  the  output,  ett),  will  consist  mostly  of  a 
sinusoid  of  frequency  f,,;  this  signal  is  then  correlated  with  the  input  to  produce  another 
sinusoid  of  the  same  frequency  at  the  summing  junction.  When  the  loop  is  closed,  the 
output  sinusoid  resulting  from  the  convolution  of  the  input  and  the  code  waveform,  s(t).  is 
1(0  out  of  phase  with  that  which  is  produced  by  the  feedback  signal,  and  hence, 
cancellation  of  the  sinusoid  takes  place.  Again,  the  feedback  signal  is  negligible  for 
small  amplitude  broadband  signals. 

Both  the  Passive  and  Active  systems  require  convolution  and  correlation  as  the  basic 
building  blocks.  in  choosing  the  architecture  to  be  used,  we  consider  the  following 
important  requirements:  l)wide  bandwidth.  2 tdynami ca 1  ly  controllable  convolver  and 
correlator.  3>convolver  -  correlator  compatibility.  4)large  dynamic  range.  These 
considerations  lead  us  to  choose  a  space-integrating  acoustooptic  architecture. 


Pig. 3  Space-Integrating  Acoustooptic  Convolver 

Convolution  can  be  performed  with  two  AOD's  as  shown  in  Fig. 3.  The  diffracted  signal 
from  the  first  AOD  is  imaged  with  unity  magnification  onto  the  second  and,  in  the  coherent 
realization  shown,  the  diffracted  components  are  spatially  integrated  onto  a  single 
detector  whose  output  photocurrent  is 

T/2 

I,(t)  •  /  ij  <  t»t ) i2 ( t-t )dx ,  t-x/v,  T-W/2.  (1) 

-T/2 

ijft)  and  ij(t)  are  the  complex  amplitudes  of  the  AOP  input  signals,  v  is  the  acoustic 
velocity  in  the  AOD.  and  K  iE  the  physical  length  of  the  AOD  aperture.  This  integral  can  be 
manipulated  to  yield  the  more  familiar  form: 


lj(t)  *  / 


i.  (t )  i.  (2t-Odx  . 


which  is  recognized  to  be  a  finite  window  convolution,  time-compressed  by  a  factor  of  two. 
The  seemingly  troublesome  compression  is  actually  an  advantage  in  this  case  as  will  be 
apparent  when  we  consider  the  correlator  implementation. 


Correlation  can  be  performed  with  the  convolution  configuration  just  described  by  tirr,e- 
reveraing  the  input  to  one  of  the  AOD'a.  This,  however,  is  not  acceptable  for  real-tin, e 
operation.  Shown  in  Fig. 4  ia  a  epace-integ rat  1 ng  processor  which  provides  the  relative 
motion  between  the  two  Input  signals  that  is  required  for  correlation  by  using  a 
demagnifying  system  between  the  two  AOD's.  The  diffracted  light  from  the  first  AOD  is 
imaged  with  a  2 : 1  magnification  onto  the  aecond  AOD.  and  the  diffracted  light  from  both 
AOD's  is  spatially  integrated  onto  a  single  detector.  If  the  input  to  the  second  AOD  is 
time-compressed  by  a  factor  of  two.  i.e..  l4(t)  ■  ij(2t>,  then  the  correlator  yields 

T/4  t*T/2 

Ij(t)  «  /  i«(t*2T>i4(t»T)dr  *  /  i j < t ) i - (t *t ) dt .  (3) 
-T/4  t-T/2 

as  its  output.  The  above  is  seen  to  be  a  finite  window  correlation  of  the  two  signals. 
ij(t)  and  ij(t>.  The  advantage  of  obtaining  the  time  compressed  output  from  the  convolver 
is  now  clear,  since  it  can  be  used  as  the  input  to  the  correlator  just  described  and  thus 
obtain  a  consistent,  cascaded  convolution-correlation  operation. 


Fig. 4  Space-Integrating  Acoustooptic  Correlator 

The  optical  implementations  of  the  Passive  and  Active  Processors  are  shown  in  Figs.  5 
and  6.  respectively,  and  they  differ  only  in  their  electrical  interconnections.  Fig.  5(f) 
is  the  system  of  Fig. 1(2)  with  the  optically  implemented  blocks  of  Figs.  3  and  4.  In  both 
cases,  the  upper  branch  of  the  processor  computes  the  correlation  while  the  lower  performs 
the  convolution.  Since  the  convolver  and  correlator  have  a  common  input,  the  first  ACD  is 
shared.  A  more  detailed  description  can  be  found  in  the  previous  publ icaticn  . 
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Fig . (  Adaptive  optical  Processor  (Active) 
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A  general  adaptive  scheme  is  shown  in  Fig. 7  in  which  each  channel  corresponding  to  one 
array  element  is  filtered  by  a  transversal  filter  with  n  controllable  tap  weights;  if  the 
array  has  N  antenna  elements,  then  the  array  processor  requires  the  adaptive  control  of  nN 
independent  weights.  Therefore,  the  computational  load  for  broadband  adaptive  phased  array 
processing  is  much  heavier  than  that  for  adaptive  temporal  filtering  because  here,  the 
array  processor  must  adapt  in  two  independent  dimensions:  space  and  time. 


Pig. 7  General  Adaptive  Array  Processing  Scbeae 


Optical  processors  are  particularly  attractive  candidates  for  adaptive  broadband  array 
signal  processing  because  of  their  inherent,  multidimensional  processing  capability. 
Moreover,  the  broadband  requirement  can  be  met  by  the  use  of  broadband  multi-channel  AOD's; 
the  large  time-bandwidth  products  that  are  available  translate  to  the  possibility  of 
implementing  very  high  order  transversal  filters.  The  acoustooptic  processor  which  we  now 
describe  is  an  extension  of  the  active  temporal  procesaor  described  earlier,  to  the  space- 
time  domain.  Bere,  we  employ  a  combination  of  multi  and  single  channel  AOD’s  to  perform 
the  required  operations. 

The  output  of  a  general,  space-time  filter,  with  a  finite  accumulation  time  can  be 

expressed  as 

N  T/2 

y(t)  "  2  ^  Xj(»)hj(t-t>  dr.  (*) 

i-l  'T/2 

where  x2<t)  is  the  signal  from  the  ith  antenna  element,  and  h^(t)  is  the  filtering  function 
in  each  channel.  A  similar  expression  is  obtained  for  the  optically  implemented  space-tire 
filter  using  two  multi-channel  AOD's  shown  in  Fig.*.  The  output  of  this  filter  1e 

N  T/2 

y(t)  -  5  (i/T )/  dt.  (i) 

i-l  -T/2 

and  we  see  that  the  only  difference  from  the  general  filter  (Eq.4)  is  in  the  tine 
compression  of  the  output.  It  can  be  shown  that  the  optimum  choice  for  h^(t)  in  Eq.(3) 
satisfies  the  system  of  integral  equations: 

N  T/2 

3  /  h! ( t*r ’ ) r < « ( t-t ' )  dt '  •  Xs ,  ( t-t ) ,  (t) 

-T / 1  J  J  1 

j-1  T/ 1 

where  Sjtt)  is  the  desired  signal  vector  and  Tjj(t)  is  the  covariance  matrix  of  the  input 
noise  giver,  by  J 

t^(t)  “  E  ln2 ( t )n j (t-t ) ) .  (2) 

ni<ti  is  the  noise  vector  appearing  at  the  array  elements. 
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Pig.*  Acoustooptic  Space-Time  Pilter 

For  adaptivity,  we  need  to  calculate  and  continuously  update  the  filter  function  h^lt) 
to  drive  the  output  to  the  optimum  result.  As  with  the  previously  described  Active 
Processor,  the  output  must  be  correlated  with  the  input  to  produce  the  appropriate  filter 
function.  Since  the  array  processor  has  n  inputs  and  one  output,  this  requires  that  we 
correlate  n  signals  with  a  common  one.  This  can  be  achieved  with  the  arrangement  shown  in 
Fig. 9  which  shows  the  use  of  a  mul ti -channel  AOD  driven  by  the  n  antenna  element  outputs, 
in  conjunction  with  a  single  channel  AOD  which  is  driven  by  the  array  output  signal. 
Specifically,  the  output  of  the  ith  element  of  the  linear  detector  array  of  the  correlator 
isgivenby 

T/2 

ti(t)  -  (l/TI  /  j(tM)x,U*2r)(Jr.  <»> 
-T/2 

For  proper  correlation  to  appear  at  each  output,  the  signal  driving  the  single  channel  AOD 
must  be  time-compres6ed  by  a  factor  of  two.  This  is  indeed  the  case  for  the  syster 
described,  and  thus,  the  AOD  implemented  space-time  filter  and  the  n-channel  correlator 
with  a  single  reference  are  compatible. 

Shown  in  Fig. 10  is  the  array  processor  system  diagram  that  shows  the  interconnections 
that  are  required,  it  is  a  direct  extension  of  the  Active  Processor  to  2-dimensions.  Thr 
output  from  each  antenna  element  is  correlated  with  the  processor  output  to  produce  the 
filter  function  for  that  element.  The  steering  vector.  Sj(-t).  determines  the  lock 


van 


direction  of  the  array  and  also  is  the  temporal  reference  signal  uaed  for  the  detection  of 
the  desired  aignal.  a(t).  It  can  be  derived  by  feeding  a  tapped  delay  line  with  s(t>;  the 
output  from  the  ith  tap  correaponds  to  a^lt)  and  the  tap  apacing  is  adjuated  to  obtain  the 
proper  look  direction: 

Sj(t)  •  s(t-ii)  •  s(t-i (d/c)coa«) ,  <9> 

where  •  is  the  desired  look  angle  from  boresight.  d  is  the  array  element  spacing,  and  c  is 
the  speed  of  light. 
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Pig. 10  Broadband  Adaptive  Array  Processor 

Fig. 11  shows  the  optically  implemented  adaptive  array  processor  with  the  AOD-impl ementec 
space-time  filter  in  the  upper  branch  and  the  n-channel  correlator  in  the  lower  one.  By 
combining  Eqs  (7)  and  <i).  the  equation  that  determines  the  filter  function,  h.(t)  is  seer, 
to  be  1 

.  T  T/J 

hi(t)»si(-t)-(C/T J)>//  »!  (t*2» )*.(t-T«» )h.  (t‘T*p )d»dx  .  (ID 

i-1  -*'*  3  3 

where  G  is  the  feedback  gain.  Under  conditions  of  low  input  SNR  and  large  feedback  gam4, 
Eq.(lO)  can  be  transformed  to  the  frequency  domain  to  yield 

N 

<  G / 4  n  3 )  J  N* (wj  )N  («, )sincJ 1 (.-.j )T/ (2» )  Jdw,  =  S*(w). 
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For  comparison,  consider  the  Fourier  Transform  of  Eq.6;  when  the  accumulation  time  is 
Infinite,  the  optimum  filter  equation  ia  given  by 

N 

3  J  *  i*sj(»).  (12) 

j-i 

rAj<«>  -  F.T.  (  Ttj<0  I. 

where  f<^(o)  is  the  spectral  density  matrix.  Identifying  the  Integral  in  Eq.li  as  the 
smoothedJestimate  of  the  spectral  density  matrix  of  the  input  noise  vector.  Fo.  11  is 
approximately  equivalent  to  Eq.12.  The  effect  of  the  finite  time  Integration  window  is 
seen  in  the  smoothing  of  the  noise  spectrum. 


Fig. 11  Adaptive  Acoustooptic  Phased  Array  Processor 
A£jUU2XledS£JI)£C£ 


This  work  is  supported  by  the  Rome  Air  Development  Center  and,  in  part,  by  the  Air  Force 
Office  of  Scientific  Research. 

iiU  Kcfei Mete 

1.  R.Riegler  and  P.  Compton .  Jr. ,  'AnAdaptive  Array  for  interference  Rejection, 

Proc.IEEE  61,748  (1*73). 

2.  J.  F.  Rhodes . 'Adapt  i  ve  Filter  with  a  Time-Domain  Implementation  Using  Correlation 

Cancellation  Loops.'  Appl.Opt.22  .  282  (1*83). 


3.  A.  Vanderlugt.  ‘Adaptive  Optical  Processor.'  Appl.0pt.2i.  «ooj  (1*82). 

4.  D.  Psaltis  and  J.  Hong,  'Adaptive  Acoustooptic  Filter.’  Appl.Opt.23.  3475  (1*84). 


APPENDIX  C 


Acomtooptic  MaptlTt  Blanil  r  tree  ulna 

John  Bong  and  Demetri  Psaltis 

California  Institute  of  Technology.  Department  of  Electrical  Engineering 

Pasadena.  California  91125 


Abatract 

Acoustooptic  methods  for  adaptive  filtering  of  temporal  signals  are  discussed.  Two 
specific  architectures  are  presented:  one  utilizing  space-integration  alone  and  the  other 
combining  both  time  and  space  integrating  techniques.  Performance  issues  regarding  the 
space-integrating  system  are  discussed  in  detail,  and  a  description  and  experimental 
results  for  the  6pace-time  integrating  filter  are  presented. 

MctgiBiiad 

Optimal  filtering  in  unknown  or  time-varying  signal  environments  requires  the  use  of 
systems  which  car.  compute  the  necessary  noise  characteristics  and  adapt  to  yield 
approximately  optimum  performance.  In  this  paper,  the  case  of  broadband  signals  received 
in  the  presence  of  strong,  additive,  narrowband  jammers  whose  spectra  are  unknown  is 
considered.  This  situation  is  illustrated  in  Fig.l.  showing  the  useful  broadband  signs! 
overwhelmed  by  a  strong  jammer.  If  the  spectral  density  of  the  jammer  were  known  apriori, 
then  the  Wiener  filter  whose  transfer  function  is  given  by1 
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would  be  used  to  yield  the  optimum  performance  in  the  least  mean  square  error  sense  (S£U) 
and  Sr(»'  are  the  spectral  densities  of  the  signal  and  uncorrelated  jammer  noise, 
respectively).  In  the  case  of  interest,  the  noise  spectrum  is  unknown  and  must  therefore 
be  computed  and  the  system  response  adjusted  accordingly  to  approximate  the  Wiener  filter. 


Fig.l  Optimum  Filtering  Fig. 2  Passive  Processor 

Several  optically  implemented  adaptive  systems  have  previously  beer,  proposed?  •  5  •  *  '  *  A 
block  diagram  of  the  adaptive  system  discussed  in  this  paper  is  shown  in  Fig. 2.  The  system 
has  been  described  in  detail  in  a  previous  publication2,  and  its  operation  is  briefly- 
reviewed  here.  The  two  major  components  of  the  system  are  the  convolution  and  correlation 
blocks  (indicated  by  •  and  *.  respectively,  in  Fig. 2).  The  outputs  of  the  convolver  and 
correlator  blocks  are  related  to  their  inputs  by 


u  (  t  )  - 

6 j  I  t-t ) 

Sj  1 1  )  dr 

wit)  - 

l: 

S j  <  t*T  > 

6 J  I  T )  dt 

(2  ' 


respectively,  where  s1(t)  denote  the  inputs  to  the  blocks.  The  processor  discriminates 
between  narrowband  signals  which  correlate  well  with  delayed  versions  of  themselves  anc 
broadband  signals  which  correlate  poorly.  Thus,  with  the  loop  open  at  the  summing  node, 
the  feedback  signal  x  ( t )  for  an  input  consisting  entirely  of  a  broadband  signal  is 
negligibly  small  because  such  an  input  correlates  poorly  with  delayed  portions  of  itself. 
This  is  contrasted  to  the  case  where  the  input  is  a  sinusoid,  in  which  case  the  convolution 
of  the  input  with  itself  delayed,  followed  by  a  correlation  with  the  input  delayed  by  the 
same  amount  yields  a  sinusoid  of  the  same  frequency  and  phase  as  the  input,  when  the  Icoy 
is  closed,  broadband  components  of  the  input  feed  through  unaffected  while  suppression  cl 
correlated  portions  takes  effect.  The  following  approximate  input-output  relaticnshii  ca: 
be  derived  for  the  filter  in  Fig.2  . 


zit)  *  xlt)  -  C//.Z  z(t)n*(«-«)n(t4o-t-o)  dad*. 
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whet*  *(t)  la  the  output.  n(t>  is  the  input  jammer  noise,  and  x(t>  •  s(t)  *  n  ( t )  is  the 
total  input.  A  aimple  interpretation  ia  found  by  the  Fourier  transformation  of  the  above 
equation,  reaulting  in  the  following  'transfer'  function: 


X (•)  "  1  ♦  G  lK(«) 1 2 

where  !(•),  X<«),  and  N(«)  are  the  Fourier  transforms  of  z(t).  *(t).  and  n(t), 

respectively.  If  the  jammer  spectrum  is  sharply  peaked  at  certain  frequencies,  the 
processor  serves  to  suppress  these  components  by  an  inverse  filtering  operation. 


Space  Integrating  IS1)  Optical  Processor 
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An  acoustooptic  convolver  can  be  implemented  using  the  SI  architecture  shown  ir.  Fig. 3. 
The  specific  set  up  shown  in  Fig. 3  is  a  coherent  version  of  a  geometry  which  can  just  as 
easily  be  implemented  using  incoherent  techniques.  The  light  diffracted  frorr  A0D1  is 
imaged  onto  AOD2 ,  resulting  in  two  counter-propagating  signals  at  the  exit  plane  of  A0D2. 
These  field  amplitudes  are  then  spatially  integrated  by  lens  L3,  and  the  integrated  signal 
i6  detected  by  the  photodetector  to  yield  the  photocurrent,  with  a  component  proportional 
to : 

T/  4 

I.(t)  •  <2/T)  i.(t*x>  i,(t-x)  dx.  T-K/v.  x-x/v,  ($) 

where  v  is  the  acoustic  velocity  in  the  AODs,  and  K  is  the  aperture  of  the  AOD.  A  change 
of  the  integration  variable  leads  to  a  more  familiar  form  given  by 


Ij(t)  -  (2/T)  /^_T/J  i x < -c )  ia(2t-x)  dx. 


which  is  recognized  to  be  a  finite  window  convolution,  time-compressed  by  a  factor  of  twc. 
The  complication  of  this  time  compressed  output  becomes  an  advantage  when  considering  the 
implementation  of  the  correlator. 
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Fig. 3  Space  Integrating  Convolver 
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Fig. <  Space  Integrating  Correlator 


Correlation  can  be  performed  by  the  set  up  shown  in  Fig. 3  if  one  of  the  input  signals 
is  t ime- r e v e r sed.  Since  this  is  incompatible  with  real-time  operation,  a  different 
geometry,  one  that  is  compatible  with  the  time-compressed  output  of  the  convolve;,  is 
examined.  Correlation  requires  that  the  input  signals  slide  past  each  other  without 
coordinate  inversion.  This  is  achieved  by  the  imp>l  ementat  l  on  shown  in  Fig. 4  ir.  which  the 
diffracted  light  frorr  the  first  AOD  is  imaged  with  a  23  demaar.  if  icatron  ente  thi  ■  ccr.c 
AOD.  At  the  exit  plane  of  the  second  AOD,  the  two  acoustic  signals  are  co-propagati  ..  i 1 1 
different  velocities  because  of  the  demagnification.  These  signals  are  then  sia:  1} 
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integrated  onto  a  single  detector.  If  the  input  signal  to  the  aecond  AOD  is  time 
compressed  by  a  factor  of  two.  i.e..  i4(t)  •  ij(Jt),  then  the  correlator  output  is 

I2(t)  •C2/T)/_^^  iJ(t*2»)i4(t*t)dT  «  (2/T)/***^*  1 3 (t > *s (t*T )dt •  (7) 

The  above  is  the  desired  finite  window  correlation  of  the  two  signals.  ij(t)  and  i3(t>. 
Since  the  output  of  the  convolver  in  Fig.)  is  Indeed  time-compressed  by  a  factor  of  two. 
it  can  be  used  to  drive  the  second  AOD  of  the  correlator,  resulting  in  the  cascaded 
correlation-convolution  operation  that  is  required  in  the  passive  processor. 
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Fig. 5  Space  Integrating  Adaptive  Optical  Processor 
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Fig.*  Freq. Response  to  One  Jammer 


Shown  in  Fig.J  is  the  complete  processor.  Since  the  correlator  and  convolver  have  one 
common  input,  one  AOD  can  be  shared  by  the  two  blocks,  resulting  in  a  system  requiring  only 
3  devices.  The  input-output  equation  describing  the  processor  in  Fig.J  is  given  by 

ztt)  •  x(t)  -  (4G/T 1 )  f  f'L/*  x*  <  t*2v )  x  ( t*x*v )  z  ( t-e«v  )dtdv .  (8 

-T/4 

where  x(t)  is  the  input  to  the  system  and  z(t)  the  output.  For  x(t)  •  s(t)  ♦  n(t).  where 
n(t)  is  the  narrowband  jammer  and  s(t)  the  broadband  signal,  the  above  equation  car  be 
Fourier  transformed  to  yield 

2<w>  i 

-  -  - - -  .  (9) 

X<»)  l  *  (G/f2)  f.z  lN(«) lJsincJl(w-«)T/4n)da 

The  above  'transfer’  function  has  the  form  of  the  inverse  filter  described  earlier,  except 
that  the  noise  spectrum  is  smoothed  by  a  sine  function  due  to  the  finite  time  windows  cf 
the  AODs . 
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The  optical  system  described  thus  far  was  implemented  for  preliminary  studs.  Jarrt-: 
suppression  was  successfully  observed;  however,  the  maximum,  static  null  depth  obtained  war 
limited  to  7dB.  It  is  important  to  theoretically  characterize  the  performance  limits  cl 
the  processor  in  te-m6  of  system  parameters  such  as  the  AOD  aperture  size  and  detector 
noise  before  we  attempt  to  systematically  improve  the  performance  level.  Three  important 


issues  tit  discussed:  1)  the  resolution  with  which  the  processor  discriminates  against 
jaamers.  2)  null  depth  limitations  due  to  detector  noise,  and  3)  transients  in  the  system 
response  due  to  the  delays  inherent  in  the  systea. 


The  resolution  of  the  system  is  characterized  entirely  by  the  size  of  the  AOD 
apertures  for  they  determine  the  length  of  the  signal  sample  at  any  one  time.  when  the 
input  consists  of  a  single  sinusoid  at  frequency  and  amplitude  A.  the  exact  form,  for  the 
Fourier  transform  of  the  output  can  be  determined  from  Eq.S  to  be 

XU) 

ZU)  .  - - - - -  ,  (10) 

1  ♦  GAJ  sine3 ( U-»0>T/4n] 

where  ZU)  and  XU)  are  the  Fourier  transforms  of  the  output  and  input  signals, 
respectively.  A  plot  of  ZU)  for  GAJ  -  I  ,T  -  10  (isecs.,  and  »0/2r  -  10  MHz  is  shown  in 
Fig.  *.  The  sidelobe  structure  is  entirely  due  to  the  finite  apertures  of  the  AODs.  The 
sidelobes  can  be  reduced  at  the  cost  of  widening  the  null  at  the  jammer  frequency  by 
introducing  apodization  functions  through  the  appropriate  illumination  of  the  AODs.  The 
issue  of  apodization  will  be  examined  more  closely  later  in  the  discussion  of  transients 
and  stability. 
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Fig. 7  system  Noise  Model 


While  the  finite  apertures  determine 
the  resolution  of  the  processor,  the 
maximum  nulling  depth  is  governed  by- 
other  mechanisms,  in  order  to  better 
understand  this  limitation.  the 
deterioration  in  the  cancellation  of  one 
jammer,  n(t)  «  Aexp(j«.0t),  due  to  the 
system  noise  contributed  by  the  detectors 
and  amplifiers  is  now  studied.  The 
system  model  depicting  the  noise  sources 
is  shown  in  Fig. 7,  where  nj(t)  is  the 
noise  from,  the  feedback  amplifier  while 
n2(t)  and  n3(t)  correspond  to  those 
arising  from  the  photodetectors  in  the 
convolver  and  correlator,  respectively. 


The  noise  processes  are  assumed  tc  te 
independent,  stationary,  complex  gaussian  processes3,  each  with  variance  a .  Since  the 
input  is  a  single  frequency  jammer,  the  convolver  and  correlator  behave  as  narrow-bandpass 
filters  centered  at  that  frequency.  If  the  space  bandwidth  products  of  the  convolver  ar.c 
correlator  are  larger  than  the  feedback  gain,  the  contribution  of  the  noise  sources 
n j  (t ) . n2 (t ) .  filtered  by  the  convolver  and  correlator  is  negligible  compared  to  nj(t>  which 
is  added  directly  onto  the  feedback  signal  x(t).  The  portions  of  noise  which  are  fed  back 
are  suppressed  by  the  correlator  and  convolver  and  can  be  neglected  to  first  order.  Thus, 
only  the  noise  due  to  ns(t)  appears  directly  at  the  output.  We  expect  that  the  performance 
limit  set  by  the  system  noise  would  not  be  reached  unless  the  jammer  amplitude  is 
suppressed  to  «,  the  noise  level.  In  our  experiments,  however,  the  jammer  in  the  output 
signal  was  clearly  stronger  than  the  system  noise,  showing  that  the  noise  performance  licit 
was  not  reached. 
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Fig. 9  Delay  Modified  Convolver  and  Correlate: 


Amore  severe  limitation  can  arise  from  the  time  delays  that  exist  in  the  AODs.  Show: 
in  Fig.i  is  an  AOD  with  the  spatial  coordinate,  x.  centered  on  the  AOD.  If  an  electrical 
signal  s(t)  is  used  to  drive  the  AOD.  then  the  amplitude  of  the  acoustic  wave  in  the  A  CD  :  s 
proport ional  to  s(t-x/v-T/2).  where  T  is  the  total  aperture  time.  This  effective;, 
introduces  time  delays  to  appear  in  the  feedback  path  which  can  cause  instabilities  in  the 
processor.  To  minimize  the  delays,  the  AODs  in  the  convolve:  and  correlate:  car.  be  shifte. 
lr  space  as  described  in  Fig. 9  to  result  in  the  following  convolve:  and  correlate:  in;„: 
output  equations,  including  the  delays 

T/  4 

convolver:  u(t)  •  U/T)  /  i,  tt-T/2*x )  l,  (t-v )  dx  .  on 
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correlator:  v(t>  *  <4/T)  i« (t-T/2^2J) i- )  dp. 

The  resulting  modified  aysteir  is  described  by  the  following  equation: 

x(t)  -  x(t>  -  1«<G/Tl>  /° ,  lTJ*  x*(t-T/i*2P)x(t-T/2*t«>)z(t-x«>)didp.  (12) 

-T/4  0 

For  a  single  jammer  input,  x(t)  •  Aexp(j.0t),  the  above  can  be  transformed  to  give  the 
following,  exact  frequency  domain  result: 

XU0> 

I(„)  .  .  - - - . -  (13) 

1  ♦  GA2  sine2 l (•-•p)T/ln)expt- j (•-■0)T/4) 

The  phase  factor  present  in  the  denominator  is  due  to  the  delays.  The  conseq  snee  of  this 
phase  factor  is  that  the  denominator  is  no  longer  positive  for  arbitrarily  large  values  of 
G.  In  fact,  the  denominator  vanishes  for  (w-w0)T/4  -  x  with  G  *  Gc  ■  nJ/4A2  .  indicating 
instability  for  gain  values  larger  than  Gc  which  corresponds  to  an  amplitude  suppression  of 
only  lOdB.  This  value  of  maximum  suppression  is  verified  by  computer  simulations  of  the 
processor  response  to  a  step  input.  Fig. 10  shows  the  envelope  of  the  response  of  the 
processor  to  a  unit  step  jammer  input,  for  G*l.  and  G-2.5,  which  is  near  the  critical 
value.  Gc.  The  initial  constant  portion  of  the  response  is  due  to  the  delays  in  the 
convolver  and  correlator.  Indeed,  the  second  plot  shows  an  oscillatory  behavior,  indicating 
marginal  stability. 


Fig.  10  Step  Response 


Fig. 11  Step  Response  with  Exponential  Apodization 

The  low  critical  gain  stems  from  the  rather  large  sidelobes  of  the  6inc2  function 
coupled  with  the  phase  factor  due  to  the  delays.  At  the  cost  of  loss  in  resolution,  the 
sidelobes  can  be  reduced  by  apodizing  the  aOD  windows,  resulting  in  systems  which  car 
tolerate  higher  feedback  gains  while  remaining  stable.  As  a  simple  example,  the  effect  of 
an  exponential  tapering  function  was  evaluated  by  computer  simulations.  Showr  ir.  Fic.ii  is 
the  step  response  obtained  with  exponential  apodization.  Convergence  is  still  rapid,  anti 
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the  suppression  is  improved  to  20dB.  Other  tspering  functions  stay  yield  stable  nulling 
performance  to  match  that  predicted  by  the  system  noife  considerations  as  discussed  earlier 
but  at  the  cost  of  loss  in  resolution. 

Space-Time  latwtlnq  Optical  PiPcciiPi 

in  the  SI  system,  the  desired  correlation  function  is  formed  by  a  SI  correlator  and 
read  out  with  a  SI  convolver.  Alternative  approaches  where  integration  in  both  tine  and 
space  are  U6ed  have  been  described  by  Rhodes9  and  Penn6.  In  such  systems,  the  correlation 
integral  is  computed  by  a  time-integrating  spatial  light  modulator  (SLM)  and  read  out  with 
a  space-integrating  convolver.  We  now  consider  an  implementation  using  the  photoref ractive 
crystal  BSO  as  the  time-integrating  SLR.  The  advantages  of  this  implementation  over  the 
previous  ones9'6  which  used  the  Hughes  liquid  crystal  light  valve  and  a  phosphor  screen, 
respectively  as  the  time-integrating  SLMs  are:  l)the  system  response  speed  is  easily 
controllable  since  the  response  time  of  the  photoref ractive  effect  depends  on  the  exposure 
level.  Da  coherent  implementation  is  possible.  3)higher  resolution  is  exhibited  by  the 
photoref  ractive  devices.  Jlsince  the  crystal  does  not  respond  to  DC  or  bias  levels,  the 
implementation  will  not  suffer  from  bias  build-up  problems  common  among  time-integrating 
systems.  4)the  photoref ractive  devices  are  simpler  since  they  require  only  the  crystal  and 
a  voltage  source  to  supply  the  external  field  which  i6  required  for  some  crystals. 

The  photoref ractive  effect  can  be  understood  in  the  following  way.  Two  writing  bears 
of  light  which  are  spatially  modulated  with  information  intersect  in  the  crystal  to  form  an 
intensity  grating  pattern.  Because  of  the  spatially  varying  intensity  pattern,  the  charge 
carriers  in  the  crystal  redistribute  themselves  into  a  space-charge  grating.  The  electric 
field  associated  with  this  charge  grating  accordingly  modulates  the  refractive  index  of  the 
crystal  through  the  electrooptic  effect,  allowing  the  read-out  of  the  written  information 
by  a  third  beam.  A  property  characteristic  of  this  effect  is  that  the  response  tire  is 
inversely  proportional  to  the  average  writing  beam  intensity  impinging  on  the  crystal. 
This  property  translates  to  the  controllability  of  the  system  response  time  in  the  present 
application,  an  added  flexibility.  In  our  experiments,  the  time  constant  was  on  the  order 
of  a  few  seconds  because  of  the  relatively  low  writing  intensities  used. 

Shown  in  Fig. 12  is  a  time  integrating  correlator  that  exploits  the  time  response 
characteristics  just  described.  A0D1  and  A0D2  serve  to  spatially  modulate  the  writing 
Argon  laser  beams  with  twc  counter-propagating  signals  Sj(t-x/v)  and  s2(t*x'v);  an  imaging 
lens  pair  is  used  to  image  the  AODs  onto  the  BSC  crystal.  The  index  grating  formed  within 
the  crystal  is  modulated  by  the  correlation  of  the  two  signals.  The  diffracted  read-out 
wave  contains  only  the  modulation  of  the  index  grating  and  is  thus  bias-free.  It  can  be 
shown  that  the  amplitude  of  the  diffracted  light  is: 


Ediff,x) 


t  <t.-t)/i 

/  e  Sj ( tj-2x /v-T/2 ) Sj ( tj*2x/v-T/2 )  dtj. 
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t-2x/v-T/2 

/  Sj  ( tj  >  62  *  t2*<  x/v  '  dt  j 

t -  2  x / v-T / 2  - t 

t  -  time  constant  of  BSO. 


which  is  the  cross-correlation  function  of  Ej(t>  and  s2(t)  with  the  spatial  coordinate,  x. 
being  the  shift  variatle.  This  correlation  can  be  detected  by  imaging  the  crystal  onto  a 
CCD  detector  array  for  serial  read-out. 
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The  incorporation  of  the  correlator  just  described  into  the  adaptive  systen  of  Fig. 2 
is  rather  simple  and  requires  the  addition  of  only  l  AOD  with  a  4:1  telescope  to  compensate 
for  the  2:1  demagnification  used  in  writing  the  correlation  onto  the  BSO.  The  system  is 
shown  schematically  in  Fig. 13.  The  output  and  input  (  z(t)  and  x(t),  respectively  ) 
constitute  the  inputs  to  the  correlator.  AOD3  is  al60  driven  by  the  input  signal  and 
diffracts  a  portion  of  the  He-Ne  read  beam  to  pass  through  the  BSO  unmodulated.  The  DC 
from  A0D3,  however,  is  Bragg  matched  to  the  index  grating  in  the  crystal  and  reads  out  the 

correlation  function.  If  all  of  the  AODs  are  driven  at  the  same  frequency,  then  the  bear 

diffracted  by  the  BSO  crystal  and  that  diffracted  by  AOD3  are  colllnear  and  interfere 

temporally  at  the  Doppler  frequency  introduced  by  A0D3.  The  two  diffracted  waves  are 

Fourier  transformed  by  lens  L J  .  and  the  resulting  intensity  is  detected  by  the 
photodetector.  The  output  current  is  the  convolution  of  the  input,  x(t),  and  the 
correlation  formed  inside  the  crystal,  and  it  becomes  the  feedback  signal,  x(t),  which  is 
then  subtracted  from  the  input  to  produce  the  output  signal.  The  resulting  input-output 
eqn.  is 

T/  g  t  * 

Z ( t )  •  X ( t )  -  G  /  /  x  <t'-2a-T/2)x(t-4o-T/2)z(t'+2a-T/2>dt'da,  (15) 

-T/8  t-t 

which  i6  similar  tc  the  SI  system  except  for  the  presence  of  both  space  and  time  integrals. 
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One  underlying  assumption  behind 
Eq.14  and  hence  15  also,  is  that  the 
crystal  is  thin  enough  to  put  the  read 
out  process  into  the  Ranan-Nath  regime. 
This  is  not  the  case,  however,  ar.c  the 
— — —  e.  input:  xit'e^O*  coupling  between  the  writing  beams  must 

'  \  be  taker,  into  account.  Coupled  wave 

*  y  y  theory  predicts  that  the  phase  of  the 

ADD  3  t — p  "  output:  index  grating  will  be  a  function  of  the 

r  a  ratio  of  the  two  write  bean  intensities, 

*\~r  t't  e-  this  was  found  to  be  the  case 

4  ,  W- •  1  i~0tr  *  ]  |  experimentally.  This  introduces  a 

V'  1  -*V - i"’  complication  in  the  present  system. 

4**-*JL,  ?  ACC  fl  "I  t  ,  because,  as  soor.  as  the  processor  begins 

.  _ ^  t0  6uPPress  an  incoming  jammer,  the 

_ /  /  1  -|  f^i  output,  z(t),  begins  to  decrease  in. 

v  \J  amplitude,  causing  the  ratio  of  the 

^£>__  ^  t  a;-'  intensities  of  the  two  write  bears  which 

x(t)e-~0  emerge  frcm  AODl  and  A0D2  tc  dinir.ish. 

This  translates  to  a  phase  change  ir.  the 
feedback  signal,  possibly  caueuc 
instability  to  appear.  This  car.  be 
remedied,  however,  by  placing  an  AGC 
Fig. 13  Space-Time  Integrating  Processor  (automatic  gain  control)  amplifier  at  the 

input  to  AOD 2.  As  adaptation  proceeds 
and  z(t)  begins  to  decrease,  the  AGC  tracks  the  drop  in  amplitude  to  keep  the  input 
of  AOD 2  at  a  fixed  value.  The  system  gust  described  has  beer,  set  up  ir.  the  laboratcr>  fc: 
preliminary  study.  Shown  in  Fig. 14  is  a  result  showing  the  spectrum  of  the  input  signal,  a 
jammer  at  70  MHz.  and  that  of  the  processor  output.  The  observed  suppression  was  ato„t  15 
dP.  Better  suppression  can  be  expected  as  better  AGC  amplifiers  are  used. 


Fig. 13  Space-Time  Integrating  Processor 


Fig. 14  Suppression  of  Jammer  (70  MHz'  with  BSO  Processor 
a)  input  spectrum  b)  output  spectrum 


Conclusion 


Two  optical  implementations  of  adaptive  signal  processors  were  described.  The  space- 
integrating  processor  was  shown  to  exhibit  very  fast  response  times  but  also  stability 
problems  which  could  be  corrected  to  a  large  extent  by  apodixatlon  techniques.  The  use  of 
such  a  processor  would  be  appropriate  in  situations  where  the  jammer  noise  is  varying 
rapidly  in  time,  such  86  in  the  case  of  blinking  jammers  or  frequency  hopped  jammers.  The 
apace-time  integrating  system  was  implemented  with  a  photoref ractlve  BSO  SLM  to  form  the 
necessary  correlation  function  in  conjunction  with  a  SI  convolver  for  adaptive  filtration. 
Although  considerably  slower  that  its  SI  counterpart,  it  has  the  flexibility  of  a 
controllable  response  time.  Future  research  in  photoref ractlve  materials  may  yield  faster 
crystals. 
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